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EXISTENCE OF SOLUTIONS FOR THE KELLER-SEGEL MODEL 
OF CHEMOTAXIS WITH MEASURES AS INITIAL DATA 

PIOTR BILER AND JACEK ZIENKIEWICZ 


Abstract. A simple proof of the existence of solutions for the two-dimensional 
Keller-Segel model with measures with all the atoms less than 87r as the initial 
data is given. This result has been obtained in [11] and [1] using different 
arguments. Moreover, we show a uniform bound for the existence time of 
solutions as well as an optimal hypercontractivity estimate. 


1. Introduction 

We consider in this paper the classical parabolic-elliptic Keller-Segel model of 
chemotaxis in two space dimensions 

(1.1) u t — Au + V • («V«) = 0, 

(1.2) A v + u = 0, 
supplemented with a nonnegative initial condition 

(1.3) u(x, 0) = uq(x) > 0. 

Here for (x , 1) £ I 2 x [0, T), the function u = u(x, t) > 0 denotes the density of the 
population of microorganisms, v = v(x , t) - the density of the chemical secreted by 
themselves that attracts them and makes them to aggregate. The system (1.1)- 
(1.2) is also used in modelling the gravitational attraction of particles in the mean 
field astrophysical models, see [2]. 

As it is well known, cf. e.g. [7], the total mass of the initial condition 

(1.4) M = J uo(x)dx 

is the critical quantity for the global-in-time existence of nonnegative solutions. 
Namely, if M < 87r, then solutions of (1.1)—(1.3) (with uq — a finite nonnegative 
measure) exist for all t > 0. For the local-in-time existence, it should be assumed 
that all the atoms of the measure u o are of mass less than 87r, see [1, Th. 2]. When 
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M > 87t, nonnegative solutions blow up in a finite time, and for radially symmetric 
solutions mass equal to 87r concentrates at the origin at the blowup time, see e.g. 
[ 6 ], 

Our goal in this note is to give an alternative proof of the local in time existence 
of solutions to (1.1)—(1.3) when uq € A4(R 2 ) is a nonnegative finite measure with 
all its atoms of mass less than 87 t. We believe that this approach is conceptually 
simpler than that in the recent paper [1] (which used elaborated arguments for 
interactions of solutions emanating from localized pieces of initial data), and those 
in previous papers [10], [11]. The latter approaches used heavily the free energy 
functional for system (1.1)—(1.2) considered in bounded domains. Moreover, our 
condition (1.5) seems to be more clear, and shows that measures with small atoms, 
which are not well separated as it was assumed in [1], are also admissible as initial 
data for the system (1.1)—(1.2). Compared to [1] here, however, we obtain neither 
the uniqueness property of solutions nor the Lipschitz property of the solution map. 

The main result of this paper is 


Theorem 1.1. Let 0 < uo € L 1 
(1.1)-(1.2) such that 


be a smooth initial density for 


(1.5) |K * ls(i)||oo < 8tt - e 0 

for some fixed £o > 0 and the unit ball B( 1) centered at the origin of in R 2 . Then, 
there exists a solution of the problem (1.1)-(1.3) on the interval [0, to] with to = 
to(eo,M) such that 

(1.6) sup ||u(t)|| p < B , 

0<t<i o 

where the constant B depends on M and £q (in particular, B does not depend on 

IKIloo)- 

Note that the condition (1.5) reads 

J u 0 (y)dy= j u 0 (x-y)dy < 8tt - £ 0 

B(x, 1) B(l) 

for all the balls B(x, 1) of radius 1 centered at an arbitrary x € R 2 and this, in 
particular, means that if, more generally, uq were a nonnegative measure then its 
atoms will be strictly less than 8n. The key property of our estimate of to is that 
it depends only on M and Eo for all uq satisfying (1.5) with a given M in (1.4). 

We recall, that for each A > 0 and each solution u of (1.1)—(1.2) of mass M the 
function 

(1.7) u\(x,t) = X 2 u(Xx, X 2 t) 
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is also a solution, with its mass again equal to M. 

Of course, by a suitable scaling (1.7) of initial data we see that we may satisfy the 
assumptions of the result on the local existence in Theorem 1.1 for any nonnegative 
uo € A4(R 2 ) with its atoms strictly less than 8tt. Then, it is clear that we arrive 
at the following corollary, cf. [1, Theorem 2]. 

Corollary 1.2. The system (1.1)-(1.2) has a local-in-time solution for each initial 
nonnegative finite measure u o with all its atoms strictly less than 8i r. 

Indeed, it is sufficient to approximate (in the sense of the weak convergence of 
measures) such a measure uq by a sequence of initial data satisfying (after the 
rescaling (1.7) with a single A > 0) all the assumptions of Theorem 1.1. This 
approximation is possible by taking, e.g., e SnA ug for any sequence S n \ 0. Then, 
the existence time to is bounded from below by a positive quantity (since depends 
on M and A only). Next, we infer from the hypercontractivity estimate (1.6) 
and from the standard regularity theory for parabolic equations that for every 
multiindex a 

\\D a u{t)\\ p < C a Bt 1/p ~ 1 ~ |a|/2 , 

which permits us to pass to the limit with (a subsequence of) the approximating 
solutions which are, in fact, smooth on R 2 x (0,to)- We obtain in such a way a 
solution to (1.1)-(1.2) with the measure uq, and this solution is also smooth on 
R 2 x (0,t o )- 

The proof of Theorem 1.1 will be a consequence of a well-known fact in [2, 4] on 
the estimate of the existence time for a solution by mass of the initial condition only, 
see 2.7 in Sec. 2, by using a rather delicate argument of localization repeatedly. 

The existence of solutions results are proved (e.g. as in [2]) for the integral 
formulation of the system (1.1)—(1.3) 

(1.8) u(t) = e tA uo + B(u, u)(t), 

whose solutions are called mild solutions of the original Cauchy problem. Here, the 
bilinear term B is defined as 

t 

(1.9) B(u,z)(t) = — j (Ve^ -S ) A ^ • (u(s) V(-A)^ 1 2 (s)) ds. 

o 

It is well known that the heat semigroup e tA , satisfies the following L q —L v estimates 

( 1 . 10 ) We^zW^Ct^-^M, 
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and 

(1.11) ||Ve tA z|| p < 

for all 1 < q < p < oo. Moreover, for each p > 1 and z G L 1 (R 2 ) the following 
relation holds 

(1.12) limt 1 - 1 / p ||e tA dL = 0. 

t- 5-0 ^ 

This is the consequence of, e.g., much more general inequality valid for every finite 
measure fi £ ,M(R 2 ) and every p > 1 

(1.13) limsupf 1_1/p ||e tA /x|| p < C'p|| / u at ||A 4 (R 2 ) = Cp |m(M)|, 

{x:([i({x}) 7^0} 

where p a t denotes the purely atomic part of the measure p. The proof of (1.13) is 
contained in [9, Lemma 4.4]. This fact, equivalent to the condition (1.5) rescaled 
to other balls of a fixed radius, see (1.14) below, is crucial in the analysis of appli¬ 
cability of the Banach contraction argument to the equation (1.8). 

We recall that the formulation of our existence results in [2] used in fact condition 

(1.13) in the definition of the functional space where solutions have been looked for: 
{u : (0, T) L?(R 2 ) : |u| = sup 0<t<T t 1 - 1 /p|| u (*)|| p < oo}, and then a smallness 
condition has been assumed on the quantity |||e tA wo|||- 

The heuristics behind the argument leading to the proof of Theorem 1.1 is the 
following: the initial data diffuse into a domain which size grows as t 1 / 2 in time as 
in Corollary 2.8. Thus, we need to find a moment of time r > 0 when a counterpart 
of the condition (1.5) 

(1.14) IKr) * l B(e) || < rn 0 

holds with a sufficiently small mo given in (2.10) and g > 0 suitably small in order 
to apply the local existence result in Theorem 2.7. 

Remark. When equation (1.2) is replaced by the nonhomogeneous heat equation 
TVt = At; + u (and thus we consider the parabolic-parabolic version of the Keller- 
Segel model), the situation seems be more complicated. For instance, if r 1, 
then there exist global-in-time solutions with M > 87t which emanate from MSo as 
(purely atomic) initial data. These are self-similar solutions which are regular and 
nonunique for sufficiently large M , cf. e.g. [3] and comments in [5]. 

We will apply in the proof of Theorem 1.1 simple (but rather subtle) techniques of 
weight functions and scalings. The core of our analysis are the uniform (with respect 
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to the initial distributions) estimates on the maximal existence time, expressed in 
terms of dispersion of the initial data. 


Notations. The integrals with no integration limits are understood as over the 
whole space R 2 : f = J. The letter C denotes various constants which may vary 

R 2 

from line to line but they are independent of solutions. The norm in L P (R 2 ) is 
denoted by || . || p . The kernel of the heat semigroup on R 2 , denoted by e tA , is given 
by G(x, t ) = (47rf) —1 exp • 


2. Proof of Theorem 1.1 

The proof of the estimate of the existence time in Theorem 1.1 is split into several 
lemmata. 

For any fixed Xo G R 2 we define the local moment of a solution u by 

(2.1) A(t) = / ip(x — xo)u(x, t) da;. 

Here the weight function 

( 2 . 2 ) 

ip(x) = (l — |aa| 2 )with Vi/’(a;) = —4a;(l — |a;| 2 ) + and Aip(x) = 16|x| 2 — 8 > — 8, 

is a fixed radial, piecewise C 2 , nonnegative function ip, supported on the unit ball 
such that ip( 0) = 1. Our particular choice of the function ip is not critical. 

Lemma 2.1. Suppose that w = w(x) is a nonnegative function locally in L 1 C\L°°, 

J w(x) da; < to, and g, 6 £ (0,1). 

B{ l) 

(i) Then there exists a number Hq € (0,1) such that J w(x)dx < (1 — S)m implies 

B(e) 

f ip(x)w(x) da; < (1 — Ho)m. 

(ii) Similarly, there exists Hi G (0,1) such that if f w(x) da; < to and f ip(x)w(x ) dx 

B( i) 

(1 — Hi)m, then f w(x) dx > (l — |) to. 
b(q) 

(in) Suppose that the inequality J ip(x)w(x) dx < (1 — H)m holds with some H G 

(0,1). Then the bound J w(x) dx < (l — ff) m holds for f3 2 < ^ 

B(0) 
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Proof The properties (i)-(iii) are simple consequences of (2.2). Indeed, 

/ ip(x)w(x) da; < / w(x) da; + sup 4>(x) • / w(x) da; 

J B(l)\B(p) J 


B(e) 


B(i)\B( e ) 


I da; 


< J w(x) da: + (1 — g 2 ) 2 J w(x)< 

B(e) B(i)\B( e ) 

= (1 — g 2 ) 2 J w(x) da; + (l — (1 — g 2 ) 2 ) J w(x)dx 

B( i) B( e ) 

< (1 — g 2 ) 2 m + (l — (1 — g 2 ) 2 ) (1 — S)m 

= (1 - H 0 )m, 


where 1 — Ho = (1 — g 2 ) 2 + (l — (1 — g 2 ) 2 ) (1 — 6) = 1 — 6 (l — (1 — g 2 ) 2 ). 

(ii) is equivalent to (i) with <5 replaced by 1 5. 

(iii) For \x\ < /?, j3 2 < -f- and H < 1 the inequality ip(x) > jffrr holds. □ 

Next, we show a result on the dispersion of the initial data evolving according 
to (1.1)—(1.2). 


Lemma 2.2. Let u be a solution to (1.1)-(1.2) such that for t G [0,A] 

(2-3) ||«o * ls(R 0 )lloo < m 

for some A > 0, Rq = 6 • 128^^ > 0 and mo < to < 8n — £q- Then, there 

exist numbers A\ = Ai{M,eo), 6 = 5(M, ecu too) and g = g(M, eo,mo) such that if 

f u(y,t)dy > (1 — S)m for some t G [0, A], then the differential inequality 
\y-x 0 \<e 

A ft) < -■& 

holds with some id = $(M, £o,?no) > 0. 


Proof First we give a uniform estimate of the time derivative of the moment Aft): 

|A'(i)| < C M . 

Let us compute the time derivative of A using equations (1.1)—(1.2) and (2.2). 
Symmetrizing the bilinear integral / u(x,t)'\7v(x,t) ■ VV’(a:) da; with the solution v 
of (1.2) given by v(x, t) = — A- J u(y, t) log \x — y\ dy we obtain 

A'(f) = ju(M) AVTOd* 

(2.4) +^~ JJ — • (x - y)u(x, t)u(y, t) da: dy. 

From (2.4) and (2.2) we immediately get A'(<) < 8 M + AM 2 . 
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Using (2.2), the bound \Vip(x) — Vtp(y )| < 4 and the relation f u(x, t ) dx < 87 r 

B(e) 

with g < 1 < 2 < i?o and , we arrive at 

(2.5) 


U<e2<-R 0 <|y| 


4^(1 ~~ M 2 ) + ~~ 4y(l — |y| 2 ) + 
|a: - y| 2 


(a; - y)u(x, t)u(y, t) dx d y 


16-7T M 

— n 4 - 
itn 


Next, in the annulus g < |y| < Rq we have 
/' /' Vt/’(x) - V^(y) 


( 2 . 6 ) 

||.r|<e e<\v\<Ro 
where we applied the bound 


\x ~ y | 2 


{x - y)u(x,t)u(y,t) dxdy 


< B8tt(5, 


(2.7) 


|V^(a;) - ViA(y)| < B\x - y| 


which holds for some constant B , as well as / u(y, t) dy < 5m < 5 8 tt. 

g<\y\<R 0 


< 647 t5 + 16p 2 877 . 


Finally, by (2.4) we have simply 

u(x, t)Aip(x) da; + 8 / u{x,t)dx 
B(e) 

Now, the crucial estimate for the bilinear integral in (2.4) is 
4a;(l — |x| 2 ) + — Ay(l — |y| 2 )+ 


U<e \v\<e 


-y I s 


• {x - y)u{x, t)u{y , t) dx dy 


l 


( 2 . 8 ) 


—4 


u(x, t dx 


\B(e) 


< B5M 2 


Here we used the following properties of the weight function 

\ip{x) — 1| = |(1 — M 2 )+ - !| = |2|a :| 2 - |x| 4 |+ < 2 |x| 2 , \Aip(x) + 8 | < 16|x| 2 , 
and an improvement of (2.7): 

IVV’(x) - V^(y) + 4(x - y)\ < Bg\x - y\ 
valid for all |x|, \y\ < g. Therefore, we get 


A'(t) < —8 J u(x, f)dx+ — 

( r \ 

/ u(x, t) dx 

B(e) 

\B(e) / 


128ttM 

Rn 


+1 6ttB5 +6477 d + 1 6g 2 8 tt . 
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Since 


u(x, t ) da; 


B(e) 


—8tt + J w(a:,f)da;j < — £q, 


B(e) 


it suffices to choose Rq = S < g 2 < therefore Rq = 6 ’ 1 ^ 7rM . □ 


6 C 2 > 


Note that a variant of Lennna 2.3 below has been obtained in [11] by a (rather 
elaborate) radial rearrangement argument of [8]. The proof we present uses only 
weight functions and localized moments defined by them. 


Lemma 2.3. Suppose that u = u(x,t) is a solution of (1.1)-(1.2) satisfying all 
the assumptions of Lemma 2.2. Then for all t > Aq = Ao(M,mo,eo) and Hi = 
H\(M,mo,£o) we have 

/ u(x, t) da: < (1 — Hi)m. 

B(e) 

Proof Let to = 0. If / u(x, To) da; < (1 — 6)m, then ti = to- Otherwise, let 
B(g) 

If \ 

Ti = inf < t < A : / u(x, r) dx = (1 — S)m > . 

{ B( e ) ) 

In order to obtain necessary estimates for ti, observe that by Lemma 2.2, for all 
t G [t 0 ,Ti] the inequality -^A(t) < — d holds. Since A(0) < m and A (t) > 0, we 
have T\ — To < m/'d. By the first part of Lennna 2.1 we arrive at the inequality 
A(tl) < (1 — H 0 )m. Next, we define either 

T 2 = inf {ti < t < A : A (t) = (1 — Hi)m} 

if this exists, or T 2 = A. Then for every t G [ti,T 2 ] we obtain A < (1 — Hi)m. 

If T 2 = A , we are done. If not, by the second part of Lennna 2.1 we infer that 

J u(x,r 2 )dx > (l — |) m > (1 — S)m. Therefore, A (t 2 ) < — d implies that 
B(e) 

A(t 2 — h) > (1 — Hi)m for a sufficiently small ft, > 0, contrarily to the definition 

of t 2 . Thus we get Aft) < (1 — Hi)m for t G [ti, A]. Consequently, by Lemma 2.1, 

we have J u(x,t)dx < (l — ) to, for t > ^ = A 0 (M,e 0 ) and /? < \H X ! 2 as 

B(0e) 

in Lennna 2.1 (iii). We denote /3g in the sequel again by g. □ 


Corollary 2.4. If u solves the system (1.1)-(1.2) on the time interval [0,AT], 
A > Aq, and satisfies the estimate 


(2.9) 


u(x, t) da; < to < 87r — £qi 


B(TV 2 ) 
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then for t £ [AqT, AT] we have f u(x, t) da; < (1 — H\)m. 

B(gT V 2 ) 

Proof This follows from Lemma 2.3 applied to the rescaled function ur(x,t) = 
Tu(T x / 2 x,Tt) which, evidently, is also a solution of (1.1)—(1.2). □ 


Corollary 2.5. If a solution u of (1.1)-(1.2) exists on the time interval [0, A] then 
for all t £ [Ao(l + Q\ + • • • + g[~ ), A) with gi = -^ the function u satisfies the 
estimate 


j «(*,0d*<(l -Hjm 


B(e i) 

as long as the inequality (1 — Hi) J m > mo holds. 


Proof It suffices to apply Corollary 2.4 to the functions u(x,t), u(x,A 0 + t), 
u(x,A 0 + q\A 0 + t), ..., rescaled with T 0 = 1, Tf = q\, T 3 = gf, T 3 = gf, , 
consecutively. □ 


Corollary 2.6. Suppose that a solution u of the system (1.1)-(1.2) exists for t £ 
[0,A] and satisfies 


i(t) * 1 


Then for any s £ [0,1] and t £ 


B( 1) ||oo 

t^A 0 ,A 
i-ei u ’ 

U(t) * 1 B (sr) 


<8n — £q. 
the inequality 

oo ^ rno , 


holds with r = g[ 


log 8-7T 

log(l-Hi) 



Proof Apply Corollary 2.5 to the rescaled solution. 


□ 


Now we recall the existence result in [2, 4] 

Theorem 2.7. There exists a (small) mo > 0 such that the condition 
(2.10) IK * 1b(i)||oo < 2m 0 

guarantees the existence of a local-in-time solution (on a time interval [0,T] with 
T = T(M)) satisfying the estimate ||u(t)|| p < Ct x ! p ~ x for p = |. 

In fact, we will use in the sequel the following immediate version of Theorem 2.7 
which takes into account the scale invariance (1.7) 
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Corollary 2.8. For each a > 0 there exists a > 0 such that the condition 

IK * 1b(<tt1/2)IIoo < 2Too 

implies the existence of the solution of (l.l)-(l.S) on the time interval [0, ar]. Here 
t > 0 is any small positive number and a can be chosen as a = ao<r 2 for some 
ao > 0 . 


Theorem 2.7 has been proved (even for sign-changing measures) in [2, Theorem 
2] (cf. also [4, proof of Theorem 2.2]) using a standard contraction argument applied 
to the formulation (1.8). 


Lemma 2.9. Suppose thatuo £ L 1 (R 2 )nL°°(K 2 ) is a smooth nonnegative function 
satisfying the condition 

IK * V’lloo < 87T - £ 0 . 


Then, the solution u with the initial condition (1.3) uq exists at least on the time 


interval 


0 , 


2 Cm 


Proof The inequality A'(<) < Cm implies that Aft) < 87 r — £q/2 for t < n = 
min 1 2 c M ’ t }’ w ^ iere T th e maximal existence time of u. By (iii) of Lemma 2.1 
we obtain the bound || u(t) * ls( ei )||oo < 87 t — ^ for all t < T\. From Corollary 
2.6 we infer that there exists (Jo = r\T^ such that ||u(t 2 ) * 1 b(< 7 0 )||oo < mo for 
< T 2 < Ti. By Theorem 2.7 the solution with uq = u(x,T 2 ) as the initial 
condition exists for t £ [ 0 , 0 x 2 ] with some a > 0 independent of X 2 . Therefore, by 
Theorem 2.7, this solution can be continued onto the interval [0,Ti + or 2 ]. This 
solution satisfies the estimate ||M(ri + ar 2 )|| p < Ct^ p 1 for each p £ [|,2]. Finally, 
a recurrence argument permits us to obtain a classical solution u = u(x,t) on the 
whole interval [0,T o ] with T 0 = T 0 (e 0 ,M), and applying once more Corollaries 2.5 
and 2.8 this satisfies the hypercontractive estimate ||u(t)|| p < for p £ [|, 2 ]. 

The extrapolation of that estimate to the whole range of p £ (l,oo) is standard, 
see e.g. [4]. □ 
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